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Abstract. Excitation functions of flow and stopping observables for the Au+Au system at energies from 
40 to 1500 MeV per nucleon are presented. The systematics were obtained by merging the results of the 
INDRA and FOPI experiments, both performed at the GSI facility. The connection to the nuclear equation 
of state is discussed. 
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PACS. 25.70.-z Low and intermediate energy heavy-ion reactions 
Projectile and target fragmentation 



25.75.Ld Collective flow - 25.70.Mn 



1 Introduction 

The study of collective flow in nucleus-nucleus collisions 
has been an intense field of research for the past twenty 
years [1,2]. At beam energies below several GeV per nu- 
cleon, it is mainly motivated by the goal to extract the 
equation of state (EoS) of nuclear matter from the quan- 
titative comparison of measurements with the results of 
microscopic transport-model calculations [3-5]. Consider- 
able progress has been made in this direction in recent 
years but the constraints on the EoS obtained so far re- 
main rather broad [5,6]. 

The results of flow measurements performed before 
1999 have been extensively reviewed in refs. [1,2]. In the 
meanwhile, a variety of new results has become avail- 
able regarding the directed [7-27] and elliptic [22-33] flow. 
These recent experiments have expanded the study of flow 
over a broader range of incident energies. New results 
became available on collective motion of produced parti- 
cles [12-15]. Several studies have focussed on balance (or 
transition) energies associated with sign changes of a flow 
parameter [20-22,28-31]. High statistics measurements al- 
lowed to explore the transverse momentum dependence of 
flow [17-19,27,28]. 

Since flow is generated by pressure gradients, it is clear 
that its quantitative study reveals aspects of the EoS. 
However, by itself, flow is not sufficient to fix the EoS. 
We need to know, as a function of beam energy, what 
density was achieved in the collision. An optimal condi- 
tion that matter be piled up to form a dense medium, is 
that the two colliding ions be stopped in the course of 
the collision, before the system starts to expand. Infor- 
mation on the stopping can be obtained by studying the 
rapidity density distributions of the ejectiles in both the 
beam direction (the original direction) and the transverse 
direction. Recently [7], the ratio of the variances of the 



transverse to the longitudinal rapidities was proposed as 
an indicator of the degree of stopping and it was found 
to correlate with flow provided the incident energy E/A 
exceeded 150 A MeV. While this flow-stopping correlation 
is only indirectly connected to a pressure-density correla- 
tion, it represents a potentially interesting constraint for 
microscopic simulations tending to extract the EoS from 
heavy ion data. 

The main purpose of this review is to present the exci- 
tation functions of flow (directed and elliptic) and of stop- 
ping in 197 Au + 197 Au collisions. This heavy, symmetric 
system has been studied with a variety of detectors in 
the intermediate energy domain throughout the last two 
decades: 



Experiment 



Reference E/A (MeV) 



PLASTIC-BALL [34-37] 150-1050 

MSU-ALADIN [38-40] 100-400 

LAND-FOPI [41] ' 400 

FOPI [28,30,51] 90-1500 

EOS [42] 250-1150 

MULTICS-MINIBALL [43,44] 35 

MSU-4tt [20] 25-60 

INDRA- ALADIN [22 , 45 , 48] 40- 1 50 

CHIMERA [46] 15 

The phase space coverage and the range of observables re- 
ported in these studies vary considerably. All these data 
sets could be and, in most cases, were indeed used for 
flow studies. However, except for the comparative study 
between the Plastic Ball and the EOS data on directed 
flow [42], and between the Plastic Ball, the FOPI and the 
INDRA data on elliptic flow [22,30], no detailed compar- 
ison has been made so far, in this energy domain, of the 
results obtained by different experimental groups with dif- 
ferent detectors. 
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In this work we will concentrate on the results obtained 
with the 47r FOPI and INDRA detector systems in experi- 
ments performed at the heavy-ion synchrotron SIS at GSI 
Darmstadt [7,18,22,28]. The covered ranges of incident 
energies were 90,4 MeV to 1.5 A GeV in the FOPI and 
AOA to 15CL4 MeV in the INDRA experiments. By com- 
bining the results obtained with the two detectors, having 
well adapted designs for the two different energy regimes, 
we were able to construct coherent systematics revealing 
a remarkable evolution of flow and stopping over a large 
range of incident energies. 

The observed agreement in the overlap region will serve 
as a measure of the absolute accuracy of the experimental 
data. We will focus on two aspects in this context, the sys- 
tematic errors associated with the unavoidable deficiencies 
of the experimental devices and on the systematic errors 
resulting from the analysis methods which are not neces- 
sarily independent of the former. Since the two detectors 
have different acceptances and the reaction mechanism 
evolves in the energy region covered by the two experi- 
ments, particular attention will be given to the problem 
of impact-parameter selection and to the corrections for 
the reaction plane dispersion, which need to be adapted 
accordingly. For the latter a new method has been devised 
and applied to the INDRA data. 

2 The detectors 

The INDRA detector is constructed as a set of 17 de- 
tection rings with azimuthal symmetry around the beam 
axis. The most forward ring (2° < 6i a b < 3°) consists of 12 
Si (300 ztm) - CsI(Tl) (15 cm long) telescopes. The angu- 
lar range 3° to 45° is covered by 8 rings of 192 telescopes 
in total, each with three detection layers: ionization cham- 
bers (5 cm of C^Fg at 50 mbar), Si-detectors (300 /zm) and 
CsI(Tl) scintillators with lengths decreasing from 13.8 cm 
to 9 cm with increasing angle. The remaining 8 rings, cov- 
ering the region 45° < 9i a b < 176°, have two detection 
layers: ionization chambers (5 cm of CzF s at 30 mbar and 
CsI(Tl) scintillators (7.6 to 5 cm). The total granularity 
is 336 detection cells covering 90% of the 47r solid angle. 

In the forward region (8i ab < 45°), ions with 5<Z<80 
are identified using the AE-E method. Over the whole an- 
gular range, isotope identification is obtained for 1<Z<4 
using the technique of pulse-shape discrimination for the 
CsI(Tl) signals. A complete technical description of the 
detector and of its electronics can be found in [47] , details 
of the calibrations performed for the GSI experiments are 
given in [48,49]. 

The FOPI detector [50, 51] is comprised of two main 
components: the forward Plastic Wall and the Central 
Drift Chamber, covering regions of laboratory polar an- 
gles of 1.2° < 9i ab < 30° and 34° < 6 lab < 145°, respec- 
tively. The Plastic Wall consists of 764 individual plastic 
scintillator units. Detected reaction products are identi- 
fied according to their atomic number, up to Z ~ 12, us- 
ing the measured time-of-flight (ToF) and specific energy 
loss. Particles detected with the Central Drift Chamber 
(Z < 3) are identified according to their mass (A) by 



using the measured magnetic rigidity and specific energy 
loss. The 3-dimensional tracking profits from a high equiv- 
alent detector granularity. At beam energies of 400A MeV 
and above, the forward drift chamber Hclitron can be em- 
ployed for mass identification of light fragments (Z < 2) 
at angles 7° < 6 lab < 29°. 

The FOPI detector has an effective granularity exceed- 
ing that of INDRA by about a factor of 4, a property 
matched to the increasing multiplicity of charged parti- 
cles with rising beam energy 1 . Both, INDRA and FOPI 
detectors are essentially blind to neutral particles, such as 
neutrons, ir° and 7's. The higher granularity is, however, 
not the only feature helping to cope with higher ener- 
gies. As the energy of the emitted particles rises, a level 
is reached where the principle of stopping the particle in 
a sensitive detecting material in order to determine its 
energy is no longer adequate because the material depth 
needed leads to a high probability of nuclear reactions 
undermining the energy measurement. To avoid this diffi- 
culty, one switches to time-of-flight and magnetic rigidity 
(in addition to energy loss) measurements: the appara- 
tus becomes larger and is no longer under vacuum. Hence 
the detection thresholds for the various ejectiles are raised. 
For the FOPI detectors this means that, e.g., at 90A MeV 
fragments with Z > 6 cannot be detected at midrapidity 
anymore. 



3 Impact parameter 

In a binary collision of massive 'objects', the transfer of en- 
ergy, momentum, angular momentum, mass etc. between 
the two partners will be strongly affected by the impact 
parameter b. As a consequence one expects to observe 
large event-to-event fluctuations due to impact parameter 
mixing. To be meaningful, a comparison of experimental 
observations among each other or with the predictions of 
theoretical simulations has to be performed for well de- 
fined and sufficiently narrow intervals of impact parame- 
ter. Generally in microscopic physics and, in particular, in 
nuclear physics, the impact parameter is not directly mea- 
surable but has to be estimated from global observablcs 
g characterizing the registered events. Global observables 
are determined using all or a significant fraction of the 
detected particles. 

The basic, so-called geometrical model assumption [52], 
underlying the association of an impact parameter b with 
an observed value g is that g changes strictly monotoni- 
cally with b allowing to postulate 

where the left (right) hand equation holds for g decreasing 
(rising) with b. The distribution da(g)/dg is determined 

1 The 47r-integrated charged particle multiplicities in central 
collisions increase from typically about 40 at 40A MeV to 95 
at 150,4 MeV and exceed 200 (with one quarter of them being 
charged pions) at 1.5 A GeV. 
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experimentally in terms of differential cross sections per 
unit of g in a minimum bias class of events, i.e. where 
a minimum number of conditions was required to trigger 
data taking. 

At intermediate to low incident energies, especially for 
E/A < 100 MeV, the literature abounds with an im- 
pressive diversity in the choice of global observables that 
have been used in attempts to select either narrowly con- 
strained impact parameters (keywords 'highly exclusive' 
or 'ultracentral') or events of special interest (keywords 
'fully equilibrated', 'fully evaporated', 'signals of phase co- 
existence'). The observables vary from very simple ones 
like proton, neutron or total charged-particle multiplicity 
to more specific ones as, e.g., participant proton multiplic- 
ity (N p ) [53,54], total (Et) [55,56] or light charged particle 
transverse kinetic energy (E^ 2 ) [57], ratio of transverse- 
to-longitudinal kinetic energy (Erat) [58,59], degree of 
isotropy of momenta (i?) [60,61], transverse momentum 
directivity (D) [62-65], longitudinal kinetic-energy frac- 
tion (E e ) [66,67], linear momentum transfer [68], total 
kinetic-energy loss (TKEL) [69,70], average parallel ve- 
locity (Va V ) [71], midrapidity charge (Z y ) [72], total charge 
of Z>2 products {Zbound) [73,74], longitudinal component 
of the quadrupole moment tensor (Q zz ) [75]. Even more 
complex observables are those obtained from sphericity 
[76,77], from the kinetic energy tensor [78,80] or mo- 
mentum tensor [67,81,82], the thrust (T) [67,83,84], the 
deflection angle of the projectile (Qdefi) [85], the flow 
angle (0fi ow ) [3,86], the location in a 'Wilczyhski plot' 
[69,86,87], harmonic moments (H2) [86,88,89], or com- 
bined global variables (p) [90]. The most sophisticated 
methods used for impact-parameter selection are based 
on, e.g., principal component analysis (PCA) [91-93] or 
on neural- network techniques (NN) [94-96] . 

There are also more technical event selection schemes 
involving the postulation of 'complete' events by demand- 
ing that nearly the full system charge or the full total 
linear momentum is accounted for. These latter methods 
are specific for a given apparatus since these observables, 
strictly constrained by conservation laws, would not be 
impact parameter selective when using a perfect detection 
system. In this case, a comparison of different experimen- 
tal data sets at a high level of precision is difficult and 
a comparison with theoretical approaches must use appa- 
ratus specific filter software that reproduces the hardware 
cuts causing the observed selectivity. In the present study, 
aiming towards joining up the data of two rather different 
setups, we will try to avoid using such concepts. We will 
restrict ourselves to the use of 'simple' global observables 
such as total charged particle multiplicity M c or trans- 
verse energy E± or its variants E^ 2 (limited to Z < 2) 
and Erat which, although it involves also the longitudinal 
kinetic energy, is highly correlated to E± due to energy 
conservation constraints. 

The quality of the achieved selectivity in impact pa- 
rameter is illustrated in fig. It shows distributions of 
the scaled impact parameter 60 — b/b max as obtained from 
the IQMD transport code [97] simulations for the reaction 
197 Au + 197 Au at 15(L4 MeV ^ We take bmax = !.i5(^V 3 + 
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Fig. 1. Simulated reduced impact parameter distributions for 
Au+Au collisions at 150^4 MeV using the global observables 
Erat (histogram) or charged particle multiplicity (crosses) for 
event selection. The two peaks correspond to nominal central- 
ities &o < 0.15 and 0.45 < bo < 0.55, respectively, as indicated 
by the vertical lines. 



1/3 

Aji ) fm and estimate b from the calculated differential 
cross sections for the Erat or multiplicity distributions, 
using the geometrical sharp-cut approximation. The fig- 
ure gives an idea of the achievable impact parameter res- 
olution, typically 1 to 2 fm for Au on Au, an unavoidable 
finite size effect. The semi-central event class, at this en- 
ergy, happens to be almost invariant against the choice 
of the selection method. For the central sample, about 
130 mb here, the Erat selection is somewhat more effec- 
tive than the multiplicity selection, an observation [30] 
found to hold for all higher energies studied with FOPI. 
We also conclude that with this selection technique cross 
section samples significantly smaller than 100 mb cannot 
be considered as representative of the chosen nominal b 
value. 

In this simulation perfect 47r acceptance was assumed. 
In reality, limitations of the apparatus will further reduce 
the achievable selectivity. For the case of FOPI, extensive 
simulations suggested that the additional loss of perfor- 
mance is small, provided the incident energy per nucleon, 
E/A, is at least 150 MeV and the considered range of 
reduced impact parameter does not significantly exceed 
6o = 0.5. 

At sufficiently high E/A, the measured directed flow 
can be used for a model-independent comparison of the 
relative performance of different selection methods. This 
is illustrated in fig. [21 with FOPI data for the reaction 
Au+Au at 400^4 MeV and for impact-parameter selections 
using either Erat or the multiplicity of charged particles 
in the geometrical sharp-cut approximation. 
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2. Mean charge- integrated (Z < 10) scaled directed flow, 
Pxdir measured with FOP1 for Au+Au collisions at 400,4 MeV 
as a function of the scaled impact parameter bo as determined 
with Erat (triangles) and with the multiplicity of charged par- 
ticles (crosses), after [7]. 



The scaled directed flow is p^ ir = p x dir/u\ cm where 
Pxdir = J2 si 9 n (y) Zu x/J2 z ( z fragment charge, u lcm 
spatial part of the center of mass projectile 4-velocity, 
u x = f3 x j is the transverse projection of the fragment 
4-velocity on the reaction plane [98]). The sum is taken 
over all measured charged particles with Z < 10, exclud- 
ing pions, and y is the cm. rapidity. For symmetry rea- 
sons, p x ° d \ r has to converge to zero as b — ► 0. The figure, 
therefore, indicates that (i) the b resolution is not perfect 
in either case and (ii) for the most central collisions the 
Erat selection provides a more stringent impact param- 
eter resolution than the multiplicity selection, as already 
expected on the basis of the simulations (fig. ^) . The max- 
imum value of p x °) ir , on the other hand, and the bo interval 
where it is located are robust observables which do not sig- 
nificantly depend on the selection method. Based on these 
observations, when FOPI data is analyzed, in general one 
employs a mixed multiplicity- Erat strategy for centrality 
selection. 

Not all global observables behave monotonically with 
impact parameter, as evident for p xc iir from fig. [3 If they 
are used to select central collisions an additional cut is re- 
quired to suppress the high bo branch. Non-monotonic be- 
haviour can also result from losses of heavy ejectiles close 
to zero degree or close to target rapidity. These losses tend 
to increase with decreasing E/A and (or) increasing bo- In 
the FOPI case we limit our analysis to bo < 0.5 and require 
that at least 50% of the total charge has been identified, 
a moderate, apparatus specific, constraint that does not 
significantly bias the topology of central collisions. 

While for particle multiplicities the idea of a mono- 
tonic b correlation is intuitively expected, this is not self- 
evident for transverse energy. At sufficiently high energy 



(> 100 A MeV), transverse energy is increasingly gener- 
ated by the repeated action of many elementary collisions 
on the nucleonic level. Since the number of such colli- 
sions increases with increasing target-projectile overlap, 
high transverse energies are correlated with low impact 
parameters. If E/A is smaller than about 100 MeV mean 
field effects involving the system as a whole dominate. 
One observes deflections of the projectile-like and target- 
like remnants to finite polar angles generating transverse 
energies that are associated with large impact parame- 
ters which carry large angular momenta. This complica- 
tion can be avoided by using the sum E 12 of transverse 
momenta of light charged particles (Z < 2) which is more 
strongly related to the dissipated energy and does not in- 
volve properties of heavier fragments. 



40 A MeV 



150 AMeV 
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Z=6 




Z=15 
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Fig. 3. Longitudinal scaled cm. rapidity density distributions 
for 40A (left) and 150,4 (right) MeV Au+Au collisions at b < 2 
fm from the INDRA experiment for selected charges as indi- 
cated. Solid histograms - multiplicity selection of the impact 
parameter, dashed - Erat, circles - E^ 2 . 



These complexities are illustrated in fig. [3] using IN- 
DRA data for Au+Au at 40 A MeV and 150 A MeV. Shown 
are charge-separated longitudinal rapidity distributions for 
central collisions, selected with three different observables, 
multiplicity, Erat and E 12 . To the extent that stronger 
yield accumulations near midrapidity indicate higher cen- 
trality, the multiplicity binning is more selective of cen- 
tral collisions than Erat at 40A MeV while at 150A MeV 
the reverse is true. This appears more pronounced for the 
cases of larger fragments shown in the lower panels. 

For the Erat and E] 2 selections in fig. [3] the centrality 
has been defined by all the relevant reaction products ex- 
cept the one of interest. This method of excluding the 'par- 
ticle of interest' (POI) from the selection criteria allows to 
avoid autocorrelations between the studied observable and 
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the one used for the estimation of centrality. On the other 
hand, the exclusion of the POI makes the observable used 
for the impact parameter selection particle dependent, i.e. 
no longer globally event dependent. This may affect the 
partitions belonging to a given centrality bin since, de- 
pending on the particle, the event may, or may not fulfill 
the criteria for a given centrality class. It has serious conse- 
quences when the autocorrelation is strong, especially for 
low energy collisions which are characterized by the pres- 
ence of intermediate and heavy mass fragments carrying 
substantial amounts of momentum. Excluding, or miss- 
ing, such a fragment unavoidably affects the measure of 
the impact parameter and increases its fluctuations. E]^ 
does not depend on the exclusion or detection of heavy 
fragments and thus is better suited for lower energies. 

On the other hand, in the case of the INDRA detector, 
the multiplicity observable does not seem to be the opti- 
mal centrality selector at high energies (fig- HO right bot- 
tom panel) where due to inefficiencies for light particles 
(multi-hits, punch throughs), this observable may admix 
less central events with higher multiplicities of fragments 
to the most central bin. Using E^ 2 as a centrality selec- 
tor avoids switching the selection method when studying 
excitation functions. As can be seen in the figure, E 12 per- 
forms similar to multiplicity at 40^4 MeV and similar to 
Erat at 150A MeV. Since molecular dynamics simulations 
confirm this observation [99] , we choose E± in the follow- 
ing as a centrality measure for the INDRA data, unless 
indicated otherwise. 



4 Rapidity density and stopping 

Rapidity distributions in longitudinal (y z ) and in an arbi- 
trarily fixed transverse direction (y x ) as obtained with the 
FOPI and INDRA detectors for central Au+Au collisions 
at 150A MeV are shown in fig. El 

To allow a closer comparison of the shapes the dis- 
tributions have been normalized to the unit area, indi- 
vidually for each fragment charge. The Erat observable 
constructed from all detected reaction products except 
the particle of interest was used as impact-parameter se- 
lector. In the case of the FOPI data, the distributions 
have been reconstructed for the uncovered phase space 
and symmetrized with respect to the cm. rapidity us- 
ing two-dimensional extrapolation methods [100] in the 
transverse-momentum vs. rapidity plane. For Z = 1,2 
these corrections represent less than 10% of the total yield, 
for heavier fragments they amount up to 30%, leading 
to estimated uncertainties of 10% near midrapidity and 
of 5% for \y\/y pr > 0.5. The INDRA distributions have 
been corrected for the 10% geometrical inefficiency [47] 
by multiplying the yields with a factor of 1.11. The po- 
sitions of the detected particles and fragments were uni- 
formly randomized within the active area of the detection 
modules. For Z = 1 the backward cm. distribution was 
used and reflected into the forward hemisphere which is 
affected by losses due to punch-through of energetic par- 
ticles. For heavier charges the forward part was used and 
symmetrized to profit from the higher granularity of the 
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Fig. 4. Yield distributions as a function of the scaled longitu- 
dinal (left) and transverse (right) rapidity for several fragment 
species within Z — 1 — 8 for central Au+Au collisions at 150^1 
MeV measured with FOPI (circles) and INDRA (histograms). 
Impact parameter bo < 0.15 is selected using Erat for both 
cases, the spectra are normalized to permit an easier compar- 
ison of their shapes. 



detector there and to avoid the higher thresholds affecting 
the yields at backward angles. 

Taking into account the systematic errors (not shown 
in the figure), the agreement of the two independent mea- 
surements is very good. This feature is far from trivial: 
due to different acceptances, especially for heavier frag- 
ments, the composition of the global event selector cannot 
be made strictly identical for the two detectors. Since at 
this incident energy the difference between rapidity and 
velocity is small, one can say that in a naive thermal 
equilibrium model, ignoring flow and partial transparency 
effects, the two kinds of distributions, longitudinal and 
transverse, ought to be equal, with the common variances 
being a measure of the (kinetic) temperature. Clearly, this 
is not the case, the transverse widths are smaller than the 
longitudinal widths, even though the selection method, 
using maximal Erat, is definitely biased towards isolating 
the event sample (on the 130 mb level) with the largest 
ratio of transverse-to-longitudinal variances (although, as 
mentioned earlier, autocorrelations were removed). 

The integration over rapidity yields absolute charged- 
particle distributions dN/dZ. The results for Au+Au at 
150A MeV are shown in fig. El For FOPI, only Z < 8 
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Fig. 5. Charge multiplicity distributions in central (fj < 0.15) 
collisions of Au+Au at 15CL4 MeV. Open circles: FOPI data; 
closed circles: INDRA data. 



yields are available at this energy while the INDRA data 
extend over almost 6 orders of magnitude up to Z = 20. 
The observed yields of heavy fragments are small, how- 
ever, only about 2-3% of the available charge is cluster- 
ized in fragments with Z > 8, as expected at this energy 
where the cm. collision energy amounts to four times the 
nuclear binding energy. With these 3% added, the FOPI 
data account for 97% of the total system charge which 
is consistent with the 47r-reconstruction method. The IN- 
DRA yields are systematically lower than FOPI by be- 
tween 10% and 30%. The lower Z — 1 yield is mainly 
responsible for the detection of only 80% of the total sys- 
tem charge with INDRA but similar differences are also 
observed for larger Z. They are most likely caused by reac- 
tion losses and edge effects in the detectors which reduce 
the effective solid-angle coverage if Z identification is re- 
quired. The light particle yields may also be affected by 
the higher multi-hit probabilities at this incident energy 
at the upper end of the INDRA regime. Extrapolating 
these observations over the full range of incident energies 
studied in this work, one may expect that reaction losses 
and the multi-hit probability are considerably reduced at 
lower incident energies for INDRA while the missing yields 
at large Z in the FOPI case will be negligible at higher 
energies for the mainly central and mid-central collisions 
that are of interest here. 

The ratio of the variances of the transverse and longi- 
tudinal rapidity distributions has recently been proposed 
as a measure of the degree of stopping reached in nuclear 
collisions [7]. The ratios obtained for central Au+Au col- 
lisions at 150A MeV, after integration over the range of 
scaled cm. rapidity -l<y<l, are shown in fig.|Hlas a func- 
tion of Z. The open circles represent the FOPI data with 
error bars which include the systematic uncertainty of the 
reconstruction procedure. The INDRA data are shown for 



Fig. 6. Ratio of transverse to longitudinal variances for central 
Au+Au collisions at 150^4 MeV as measured with the FOPI 
and INDRA detectors (open and filled symbols, respectively). 
The variances were obtained for scaled cm. rapidities in the 
range -l<y<l. The chosen selections of centrality are indicated 
in the legend. 



Z < 20 and for four different impact parameter selections 
as indicated in the figure. 

The largest ratios 0.8 to 0.9 are observed for light 
charged particles (Z < 2). With increasing fragment Z, 
the ratios decrease continuously to values of < 0.1 near 
Z = 20. In the common range of fragment Z and for the 
same impact-parameter selection (Erat, b < 2 fm), the ra- 
tios measured with FOPI and INDRA are in good agree- 
ment. The selection with Erat and b < 3 f m yields slightly 
smaller ratios as expected which, however, are similar to 
those obtained with E]^. Large transverse momenta of 
light charged particles and of fragments are apparently 
correlated. Autocorrelations are not present here because 
the particle of interest is removed from the impact param- 
eter selector (see previous section). The smallest ratios of 
variances are obtained for selections according to multi- 
plicity. 

The trends as observed as a function of Z suggest that 
the heavier fragments, even in rather central collisions, ex- 
perience less stopping than lighter ones and keep a strong 
memory of the entrance channel motion. Their transverse 
momenta seem to be, nevertheless, generated in collisions 
involving nucleons or light clusters as evident from the 
correlation with E]^. The momenta of struck nucleons 
absorbed in a cluster or the recoil momenta of nucleons 
knocked out from a cluster both contribute to their final 
momenta. Their relative weight will be smaller in larger 
fragments, consistent with the observed Z dependence. 
Overall, these observations are clearly in contradiction to 
the assumption of global equilibrium including the kinetic 
degrees of freedom. Qualitatively they agree with the pre- 
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Fig. 7. Excitation function of the degree of stopping, vartl, in 
central Au+Au collisions (b < 2 fm) obtained from the FOPI 
(open circles) and INDRA (dots) measurements. The result at 
15 A MeV corresponds to a less central selection (b < 5 fm). 



dictions of quantum molecular dynamics calculations for 
fragment production in this energy range [101]. 

A global observable to describe stopping, vartl, has 
been introduced in ref. [7] . It is defined as the ratio of the 
transverse over longitudinal variances of the summed and 
weighted rapidity distributions. The excitation function 
of this observable is presented in fig. [7| for central Au+Au 
collisions (b < 2 fm) and for the full energy range covered 
in the FOPI and INDRA experiments. 

The FOPI results have been obtained with the Erat 
selection. The data have recently been reanalyzed by tak- 
ing additional small corrections due to energy losses in 
structural parts of the detector into account. At the low- 
est three energies, this has led to an increase of vartl by 
up to about 10% compared to the data published in [7], 
while for the other energies the results are unaffected. The 
measured range of fragments extends up to Z — 6, 8, 8 for 
E = 90A, 120A, 150A MeV, respectively The contribu- 
tion of heavier fragments to the vartl observable has been 
estimated by extrapolating their weights and variance ra- 
tios to higher Z . At 90 A MeV this correction amounts to 
about 8%. The errors given in the figure are systematic 
and mainly reflect the uncertainty of the reconstruction 
procedure. 

For the INDRA central event samples, the light-particle 
transverse energy E^ 2 has been used to select b < 2 fm for 
40A to 150A MeV and b < 5 fm for the data sample at 15 A 
MeV with low statistics, originally only taken for calibra- 
tion purposes. For the charge- weighted average, fragments 
up to Z = 60 have been included. Except for the result 
at 15 A MeV, the error bars correspond to the variation 
of vartl for centrality selections within b < 1.5 fm (upper 
end of the error bar) and b < 2.5 fm (lower end) and thus 



represent the systematic uncertainty associated with the 
impact parameter determination. The statistical errors are 
below the percent level, except at 15A MeV, where they 
are the main contribution to the error shown in the figure. 

The obtained excitation function of stopping is char- 
acterized by a broad plateau extending from about 200A 
to 800A MeV with fairly rapid drops above and below. 
The highest value reached by vartl is about 0.9. With 
the INDRA data, the reduction of stopping at lower inci- 
dent energies is followed down to 40A MeV. In the over- 
lap region, a very satisfactory agreement within errors is 
observed. The measurement at 15A MeV suggests that 
stopping goes through a minimum at or below 40 A MeV. 

It is clear that only a dynamical theory will be able to 
reproduce this excitation function. Using the relativistic 
Boltzmann-Uehling-Uhlenbeck (RBUU) transport model, 
an analysis of the combined FOPI stopping and flow (see 
later) data was recently presented [102] . The input to com- 
puter codes implementing transport theoretical models are 
the nuclear mean field U (or EoS) and the nucleon-nuclcon 
(nn) cross sections a nn . Although both are not indepen- 
dent in a consistent theory, it is useful to consider their 
effects separately. In general one finds that the cross sec- 
tional part is dominant over the mean field part for a 
quantitative account of the observed incomplete stopping: 
note that if global equilibrium, or even local equilibrium 
(ideal hydrodynamics) , were valid cross sections would be 
irrelevant. Starting at the low energy end one qualitatively 
expects, when raising the energy, that the increasingly re- 
pulsive mean field (due to increasing compression) and the 
drop in Pauli-blocking of final and intermediate states in 
nn scattering (due to the increasing initial rapidity gap) 
conspire to raise rapidly the generation of transverse en- 
ergy at the expense of the longitudinal energy. At the 
higher energy end (say beyond \A GeV) again both as- 
pects (mean field and collisions) more or less may add up 
to make the drop faster. At l.hA GeV roughly one quarter 
of the nucleons are excited to a resonant state. The open- 
ing up of nucleonic degrees of freedom may lead to a soft- 
ening of the EoS. On the other hand the in-medium Dirac 
masses M£> are predicted to drop substantially in covari- 
ant theories [103-105], a fact that will seriously modify 
the phase-space and kinematical factors influencing the 
elementary cross sections [106-108]. The calculations of 
ref. [102] suggest that these in-medium modifications of 
<j nn are indeed necessary to reproduce the observed stop- 
ping. 

Besides the 'global' information shown in fig. the 
'particle differential' information reveals additional infor- 
mation on the stopping mechanisms. Figure shows that 
the partial transparency is predominantly experienced by 
the heavier fragments, which presumably have survived 
because their constituent nucleons have suffered a less vi- 
olent average collision history. This feature is also observed 
at the high energy end, although the 'heavy fragment' role 
is played there by mass A — 2 — 4 ejectiles [109]. Restrict- 
ing the stopping observable to the lightest species at the 
various incident energies, one obtains higher vartl values 
and flatter excitation functions. The combined role of the 
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mean field and of in-medium modified cross sections will 
be picked up again in section0where the flow information 
will be added to the analysis. 



5 Flow, reaction plane and corrections 

Originally, the directed flow has been quantified by mea- 
suring the in-plane component of the transverse momen- 
tum [98] and the elliptic flow by parametrizing the az- 
imuthal asymmetries using the Fourier expansion fits [110, 
111]. More recently, it has been proposed [112] to express 
both, directed and elliptic flow in terms of the Fourier co- 
efficients and i>2, respectively) and also to investigate 
the higher flow components. The coefficients v n are ob- 
tained by means of the Fourier decomposition [112-114] 
of the azimuthal distributions measured with respect to 
the true reaction plane: 

dN N ( „ v-^ , , , A 

^^)=^[ 1 + 2 ^ C0Sn ^-^) (2) 

with <j)R being the azimuth of the latter. In general, the co- 
efficients v n = (cos n {4> — 4>r)) may depend on the particle 
type, rapidity y and the transverse momentum px- 

The standard methods of measuring flow can be split 
into those using explicitly the concept of the reaction plane 
[98,112-114] and those based on the two-particle azimuthal 
correlations [115]. Still other methods have been proposed 
recently, satisfying the needs of high energy experiments: 
the 'cumulant' methods [116-118] using multi-particle cor- 
relations and the method based on the Lee- Yang theory of 
phase transitions [119,120]. The latter is expected to per- 
form well above about 100A MeV [119], while the three- 
particle variant of the 'cumulant' method is claimed to be 
useful for extracting v\ coefficients at energies near the 
balance energy and in the ultrarelativistic regime [118]. 
However, because the correlation methods require high 
event multiplicities and high-statistics data, and because 
the correlation between a particle and the flow vector 
is usually much stronger than that between two parti- 
cles [121], the reaction plane methods are still more com- 
monly used at intermediate energies. They have also been 
applied in the present case. 

Since detectors do not allow to measure the angular 
momenta and spins of the reaction products, the orien- 
tation of the reaction plane can only be estimated using 
the momenta. The resulting azimuthal angle, <pE, has a 
finite precision, and the measured coefficients v^ leas are 
thus biased. They are related to the true ones through the 
following expression [113]: 

v™ eas = (cosn(0 -4> E ))=v n (cos nA<f>) (3) 

where the average cosine of the azimuthal angle between 
the true and the estimated planes, (cos nAcfi) = (cos n(0ft— 
4>e)}, is the required correction (also referred to as 'event 
plane resolution' or just 'resolution') for a given harmonic. 



Note that, since the true values of flow are obtained by di- 
viding by the average cosine, they are always larger than 
the measured ones. 

The literature offers many different methods to esti- 
mate the reaction plane, like the flow-tensor method [78], 
the fission-fragment plane [122], the flow Q- vector method 
[98], the transverse momentum tensor [123] (also called 
'azimuthal correlation' [124]) method or others [125]. 

Among them, the Q-vector method has received spe- 
cial attention. Originally, the Q-vector has been defined 
as a weighted sum of the transverse momenta of the mea- 
sured N reaction products [98] : 

N 

Q = Y,"iPi ( 4 ) 

i=l 

with the weights lj chosen to be +(-)l for reaction prod- 
ucts in the forward (backward) cm. hemisphere and with 
the possibility to exclude the midrapidity zone. The choice 
of the optimal weights is discussed in [72,114,117,126,127]. 
Definition Q can be extended to Q-vectors built from 
higher harmonics [114], thus e.g. allowing to profit from 
strong elliptic flow, when applicable. Usually in the flow 
studies, the POI is excluded from the sum in ® to avoid 
autocorrelations. This does not concern the corrections, 
since the sub-events (see below) do not share particles. 

The corrections for the reaction-plane dispersion can 
be obtained using various methods [98,112-114,121,123, 
128-131]. What they all have in common, is the underlying 
assumption of the applicability of the central-limit theo- 
rem. In most of these methods the correction is searched 
for using the sub-event method [98], which consists in 
splitting randomly each event into two equal-multiplicity 
sub-events and getting the correction from the distribu- 
tion of the relative azimuthal angle, Z\<?i2 , between their 
individual Q-vectors ('sub-Q-vectors'). This is done either 
by using the small angle expansion [98] or by fitting with a 
theoretical distribution [113]. Instead of fitting the angu- 
lar distributions one can alternatively fit the distributions 
of the magnitude of the total Q-vector itself [112, 114]. 

Assuming the Gaussian limit, ref. [113] gives an ana- 
lytical formula for the distribution of AQn for the case 
that the distributions of sub-Q-vectors are independent 
and isotropic around their mean values. In refs. [112, 114] 
one can find the formulae relevant for the distributions of 
the magnitude of the Q-vector. 

These methods proved their usefulness for correcting 
measured flow values at higher energies (see e.g. [18,28,98, 
132, 133]) which fulfill the high multiplicity requirement. 
They are, however, not adequate for the intermediate en- 
ergy reactions, below about 100^4 MeV, where the particle 
multiplicities are lower and the events are characterized by 
a broad range of masses of the reaction products. Here, the 
applicability of the central limit theorem for devising the 
corrections is less obvious. 

Figure [8] illustrates the difficulties one encounters at 
intermediate energies. It shows the experimental distri- 
butions of A(fr\i as measured with the INDRA detector 
for the Au+Au reaction at 40 A (top) and 150 A (bot- 
tom) MeV and for two intermediate centrality bins. The 
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lines represent the fits obtained with the method described 
briefly below. The standard method [113] can, in princi- 
ple, be used to derive the corrections for energies down to 
about 80A MeV, however it fails to describe distributions 
like those at 40A MeV with double maxima or maxima 
at backward angles, which reflect the presence and impor- 
tance of the in-plane enhancement and of the correlation 
between sub-events. 
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Fig. 8. Distributions of relative angle between reaction planes 
for two random sub-events (INDRA data, dots) and fits (lines) 
using the integrated product of bivariate normal distributions 
(Eq.[SJ for centralities b ~ 5.5 — 7.5 fm (left) and b ~ 2 — 5.5 fm 
(right) and incident energies 40^4 MeV (top panels) and 15CL4 
MeV (bottom panels). 

Since at low and intermediate energies the sub-events 
are expected to be strongly correlated [129] and the dis- 
tributions of the Q-vector no longer necessarily isotropic 
[113], we have extended the method of Ollitrault [113] by 
explicitly taking into account these two effects in the theo- 
retical distribution of the sub-Q-vectors. The new method 
relies on the assumption of the Gaussian distribution of 
the flow sub-Q-vectors. This assumption has been veri- 
fied to hold even at A0A MeV, except for very periph- 
eral collisions, by performing tests with the CHIMERA- 
QMD model in which angular momentum is strictly con- 
served [134]. 

The form of the joint probability distribution of the 
random sub-Q-vectors has been searched for following the 
method outlined in appendix A of [121], by imposing the 
constraint of momentum conservation on the A^-particle 
transverse momentum distribution and using the saddle- 
point approximation. 

The resulting distribution has the form of a product 
of two bivariate Gaussians: 

d 4 N _ 1 
dQ^Q~ 2 ~ 7t%L<(i-P 2 ) ' CXP . 

(Qu - Qs) 2 + (Q 2a - Qs? - 2 P (gig - g s )(Q 2 . - Qs) 

aL(l-p 2 ) 

Qly + Q\y ~ 2 P Q\yQly 1 , . 



where we followed the convention of [113] of including the 
\/2 in a; the subscripts 1, 2 refer individually and s gen- 
erally to sub-events; the subscripts x and y refer to the 
in- and out-of-plane direction, respectively. This distribu- 
tion differs from those proposed in [113, 121, 129] in that 
it combines all three effects that influence the reaction 
plane dispersion at intermediate energies, namely the di- 
rected flow (through the mean in-plane component Q s or 
the resolution parameter \s = Qs/&sx [113]), the elliptic 
flow (through the ratio a = tJsx/<Jsy) and the correlation 
between the sub-events [129] (through the correlation co- 
efficient p e [—1,1]). It reduces to the one of [113] for 
a = 1 and p = 0, and to the one of [121] for a = 1. In de- 
riving eq. ||I>J it was assumed that the in- and out-of-planc 
correlation coefficients are equal. 

Making the division into sub-events random ensures 
that the distributions of the sub-Q-vectors are equivalent, 
in particular they have the same mean values and vari- 
ances. Since the total-Q-vector is the sum of the sub-Q- 
vectors, Q — Qi + Q2, one finds the following relation 
between the resolution parameter obtained from the dis- 
tribution of the Q-vector, and that obtained from the 
distribution of sub-Q-vectors, \s'- 

X = W2/(1 + P) (6) 

Relation ijfjj shows how the correlation between sub-events 
influences the reaction plane resolution. In particular, it 
indicates that the resolution improves in case the sub- 
events are anti-correlated (p < 0), which is predicted to be 
the case below about 150A MeV except for very peripheral 
collisions [134]. 

As in [113], the joint probability distribution (0 is 
used after integrating it over the magnitudes of the sub- 
Q-vectors and one angle, leaving the Z\<?i2 as the only 
independent variable. Unlike in [113], the resulting dis- 
tribution can not be presented in an analytical form. It 
depends on 3 parameters (x s , ct, p) which can be obtained 
from fits to the experimental or model data. The quality 
of the obtained fits is very good, even in the non-standard 
cases encountered at low energies (fig. EJ. 

The corrections for the n-th harmonic v n , depending 
now on x an d a i can be calculated (also numerically) as 
the mean values of the cos nA<f> obtained over the total- 
Q-vector distribution, in a similar way as in [113]. Figure 
shows how the elongation of the Gaussian (a) , result- 
ing from elliptic flow, modifies the corrections for the first 
two harmonics. It demonstrates that the in-plane emis- 
sions (a > 1) enhance slightly the resolution for V\ and 
considerably for V2 — even in the absence of the directed 
flow. On the other hand, squeeze-out (a < 1) deteriorates 
the resolution. The figure, in particular, shows that the 
correction can change the sign of elliptic flow in the case 
of small directed flow and squeeze-out. 

Since the correlation between the sub-events increases 
at the lower energies, the knowledge of the correlation 
coefficient p becomes crucial. Estimated values, obtained 
from model calculations, can be useful as constraints for 
p in the fitting procedure. For example, the CHIMERA- 
QMD calculations predict p to be around -0.43 for 40A 
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Fig. 9. Corrections for the first (3 upper curves) and the sec- 
ond (3 lower curves) harmonic as a function of the resolution 
parameter, x> f° r different aspect ratios, 0.9 < a < 1.5, cover- 
ing approximately the range of its variation. 



MeV and 2 < b < 8 fm and about -0.2 at 150A MeV. 
Alternatively, mean values of some rotational invariants 
which are derived from the measured data can be used to 
reduce the number of fit parameters and to constrain the 
fitting routine to search for a conditional minimum [135]. 

Instead of fitting the azimuthal distributions one can 
express the probability distribution JSJ in terms of the 
components of one of the sub-Q-vectors in the reference 
frame of the other, or in terms of the absolute values of the 
sum and of the difference of the sub-Q-vectors. The cor- 
responding 2-dimensional experimental distributions can 
then be fit using such formulae. The method of fitting the 
distributions of components of the sub-Q-vector has been 
found sensitive enough to perform well without additional 
constraints. 

The corrections obtained using various methods are 
presented in fig. They are close to one, independent 
of the method, for the range of higher incident energies 
(E > 100A MeV) where the directed flow is large and the 
reaction plane well defined by the high-multiplicity distri- 
bution of detected particles. At around 50 A MeV, they go 
through a minimum and depend strongly on the chosen 
method. The FOPI flow results, as published in refs. [18, 
28], have been corrected using the standard method, ex- 
cluding the midrapidity region of ±0.3 of the scaled cm. 
rapidity from the Q-vector to improve the resolution. The 
corrections used here for the INDRA data are obtained 
with the new method in two ways, by fitting the azimuthal 
distributions and by fitting the distributions of compo- 
nents of the sub-Q-vectors. The mean values are given in 
fig. 1101 (full circles) with error bars representing the sys- 
tematic uncertainty as given by the difference of these re- 
sults. At 15A MeV, the statistical errors dominate. Even 
at their minima, the corrections are not smaller than 0.6 
and 0.5 for directed and elliptic flow, respectively, indi- 
cating that the measured flow values will increase, after 
applying the corrections, by no more than about a factor 
of two. 
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Fig. 10. Corrections for the first harmonic and impact pa- 
rameter 6 = 2 — 5.5 fm (top) and for the second harmonic 
and b — 5.5 — 7.5 fm (bottom) used in the flow analysis of the 
Au+Au reaction. The corrections for the INDRA data (full 
symbols) are shown as obtained with the new method (dots, 
see text), with the standard (std) method (ref. [113], trian- 
gles) , and with the standard method and the restriction Z < 10 
(stars). Open circles represent the results for the FOPI data 
obtained with the standard method. 



For a comparison of the different methods and of their 
applicability, also the corrections according to the stan- 
dard method of Ollitrault [113] have been determined. 
This corresponds to fixing the parameters a = 1 and 
p = in the new method. Near 50^4 MeV, the results 
are close to zero which would require nearly infinitely 
large corrections (triangles in fig. The figure, fur- 
thermore, shows the same corrections according to the 
standard method as obtained for the FOPI case (circles). 
They are very similar and, in the overlap region, virtu- 
ally identical to the result for INDRA if the limit Z < 10 
of the FOPI acceptance is applied in the INDRA case 
(stars). This very close agreement is not unexpected be- 
cause good agreement was already observed for the un- 
corrected flow data obtained with the two detection sys- 
tems [22,28,33]. The standard method, nevertheless, fails 
below about 80^4 MeV. As mentioned above, the inde- 
pendent, isotropic Gaussian approximation is no longer 
confirmed by satisfactory fits of the experimental distri- 
butions. 
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Several additional observations and comments can be 
made. Comparing the results of the new and standard 
methods (filled circles and triangles) shows a dramatic 
improvement of the resolution obtained by taking the ef- 
fects of the correlation between the sub-events and of the 
in-plane enhancement into account. It should be stressed, 
that both these effects are responsible for the finite correc- 
tion for the directed flow at AOA MeV, near the expected 
balance energy [20]. Non- vanishing resolution, suggested 
by the new method, indicates that even here the reaction 
plane can be defined, and apparently, questions the oc- 
currence of the 'global' balance, which otherwise would 
manifest itself with the vanishing of (cosZ\</>). However, 
the finiteness of the corrections around 40A MeV may 
also partially result from the incompleteness of the exper- 
imental data and from the mixing of events with different 
centralities, which may add up to mask the signal of the 
'global' balance. The fitting procedure yields relatively ac- 
curate results for the corrections for the first two harmon- 
ics in case of the complete results of the simulations (e.g. 
2-5% accuracy for AOA MeV and 0.2-0.4% for 150A MeV 
and 4 < b < 8 fm [134]), and in particular, is able to re- 
veal the signal of the 'global' balance, but in the case of 
experimental data it will certainly return some effective 
corrections biased by the experimental uncertainties and 
inefficiencies. The effects of the latter may not necessarily 
drop out by applying eq. J2J to correct the measured ob- 
servables, but may require additional corrections. At the 
higher energies, the results of the standard and the new 
methods approach each other but, in the overlap region 
of the FOPI and INDRA experiments, the differences are 
still significant, and need further investigation. 

Comparing the less and the more complete data sets 
(stars and triangles, respectively) shows that the resolu- 
tion improves with the completeness of the data. Triangles 
represent INDRA events with at least 50% of the system 
charge collected to which an additional single fragment 
carrying the missing momentum and charge was added. 
This artificial completion of events was found important 
for peripheral collisions where, due to the energy thresh- 
olds, the heavy target-like fragment is always lost. The 
distributions of the relative angle between sub-events be- 
come then narrowly peaked at small relative angles which 
improves the resolution of the reaction plane. 

However, it is not only the reaction plane correction 
that relies on the completeness of the measured data. Also 
the measured y™ eas parameters are affected by the non- 
isotropic loss of particles due to multi-hits (INDRA) or 
unresolved tracks in high-track-density regions (FOPI) . A 
rough estimate of the correction [134] due to multi-hit 
losses for can be obtained, for segmented detectors like 
INDRA, by using the unfolded 'true' in- and out-of-plane 
multiplicities and calculating the true and measured mean 
Vi by integrating the azimuthal distribution J5J) over the 
in- and out-of-plane quadrants. The 'true' multiplicities 
can be estimated using the calculated (e.g. in a way sim- 
ilar to that of ref. [136]) or simulated (using the detector 
filter and the model data) multiplicity response function 
specific for a given detector. An analogous procedure can 



be applied also for vi\ however, due to the lack of the 
forward-backward center of mass symmetry of the detec- 
tor, the results may be less accurate. The flow parame- 
ters obtained from the INDRA data presented in the next 
section have been additionally corrected for the multi-hit 
losses using the above procedure. For v%, these additional 
corrections vary from about 7% at AOA MeV to about 33% 
at higher energies for Z = 1 and do not exceed 15% for 
Z = 2. For V2 and Z = 1 they increase from about 18% 
at 40A MeV to about 36% at 100A MeV and about 70% 
at 150A MeV, for the centrality bins in question. Within 
this simple procedure, the corrections depend essentially 
on the average of the in- and out-of-plane multiplicities 
and only weakly on their difference, that is why the cor- 
rections basically increase with the increasing multiplicity 
(thus with the centrality and incident energy). This ex- 
plains the large correction factor at 150^4 MeV. Neverthe- 
less, since v-i is small at this energy, the absolute change 
of the measured value due to the correction is small com- 
pared to that at lower energies. 

Generally, one may remark that, at energies below 
about 100A MeV, devising the corrections becomes a deli- 
cate task. The corrections are no longer those in the usual 
sense, say, of a few percent. Depending on the method, 
they may change the measured results by a large factor, 
mainly because of the smallness of directed flow around 
40v4 — 50^4 MeV. The accuracy relies in addition on the 
completeness of the data. Flow data free of reaction plane 
dispersions are, nevertheless, very desirable since they al- 
low to compare the results obtained with different detec- 
tors. They are also of great interest from the theoretical 
point of view, by permitting the direct comparison with 
the model predictions. In problematic cases, however, de- 
tailed filtering of the model results and treating them with 
the experimental type of analysis may still be necessary, if 
not for the direct comparison on the level of uncorrected 
observables, then for the reliable estimate of the system- 
atic uncertainties associated with the correction scheme. 

The effects of momentum conservation, distortions due 
to removal of the particle of interest (expected to be im- 
portant at low multiplicities (energies)) and possible cor- 
rections to the reaction plane resolution due to the detec- 
tor inefficiencies (missing part of the Q-vector) remain a 
subject for future study. 

6 Directed and elliptic flow 

The rapidity dependence of the slope of the directed-flow 
dvi/dy at midrapidity for Z = 1 and 2 particles, inte- 
grated over transverse momentum, is shown in fig. 1111 The 
INDRA data is combined with the FOPI data (published 
for Z = 2 in [18]), both measured for mid-central collisions 
with impact parameters of 2-5.5 fm and shown after cor- 
recting for the reaction plane dispersion. The FOPI data 
has been corrected using the method of [113] while the 
INDRA data has been corrected using the method out- 
lined in sect. In both data sets the reaction plane has 
been reconstructed using the Q-vector method with the 
weights lo = sign(y cm ) , excluding the POL In case of the 
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FOPI data the midrapidity region of ±0.3 of the scaled 
rapidity has been excluded from the Q-vector to improve 
the resolution. The INDRA data has been corrected for 
the effects of momentum conservation [137]. In both cases 
linear fits have been performed in the range of ±0.4 of the 
scaled cm. rapidity, except for the 15A MeV data where 
the range of ±0.55 was used. 
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Fig. 11. Slopes of directed flow dvi/dy for Z = 1 (top) and 
Z = 2 (bottom) particles integrated over px for mid-central 
collisions (2-5.5 fm). The open and filled symbols represent 
the FOPI [18] and the INDRA data, respectively. The uncer- 
tainty at 15^4 MeV is mainly statistical. The INDRA point, 
in brackets, at 150^4 MeV in the top panel is biased due to 
experimental inefficiencies for Z = 1 at this energy. 



The excitation function of the slope of the directed 
flow at midrapidity for Z = 1 changes sign around 80A 
MeV. The apparent minimum around 40^4 MeV is mostly 
suggested by the 15A MeV data point and should be con- 
firmed by other measurements. The FOPI data has been 
additionally corrected for the effects of unresolved tracks 
in the in-plane high track density region. This correction 
influences also the slope of the V\ rapidity distribution, in- 
creasing it by up to 15% for Z = 1 and up to 5% for Z = 2. 
The INDRA results have been corrected for the effect of 
multi-hit losses (see sect. |SJ). The still apparent discrep- 
ancy between the INDRA and FOPI results at 150 A MeV 
can be partially attributed to the losses of Z = 1 particles 



due to punch-through effects in the INDRA detector at 
high energies. Up to 10% of the difference may also come 
from the different methods used for correcting the reaction 
plane dispersions (see fig. 1101) . 

For Z = 2, the slope of v% is seen to rise monotonically 
with energy over the full range of 15 to 400 MeV per nu- 
cleon which is covered by the two experiments. Here, the 
agreement in the overlap region is slightly better reflecting 
the better efficiency of the INDRA detector for Z = 2 par- 
ticles. The trends observed for the uncorrected data [22] 
for v\ are preserved. Unlike in ref. [20], the excitation func- 
tion does not show a clean signature of a minimum (see 
ref. [22] for discussion). It changes sign between 50 and 60 
MeV per nucleon, in agreement with the extrapolated val- 
ues of the balance energy, Ei, a i, obtained from the higher 
energy measurements [42,138,139]. 

Negative flow is observed not only for Z = 1,2 (fig. Illjl 
but with even larger slopes also for other light fragments. 
This intriguing phenomenon has already been reported for 
the lighter systems 40 Ar + 58 Ni, 58 Ni + 58 Ni, and 129 Xe 
+ nat Sn, provided the T-plane-per-particle' method was 
used for estimation of the reaction plane [21]. For these 
systems, a balance energy has been determined by associ- 
ating it with the minima of the approximately parabolic 
excitation functions of the flow parameter which, in the 
cases of 40 Ar + 58 Ni and 58 Ni + 58 Ni, appeared at negative 
flow values. Negative flow values of light reaction products 
can indeed be measured experimentally, provided the de- 
tector is able to measure 'quasi-complete' events, including 
the heavy fragments. Then, the observed anti-flow of light 
products is measurable relative to the reaction plane fixed 
and oriented by the heavy remnants. 

A possible scenario of the anti-flow has been proposed 
for the lighter systems in [21], and for the heavy systems, 
emphasizing the role of the strong Coulomb field, in [22]. 
Despite the appeal of a globally defined balance energy, it 
is worth noticing that directed flow apparently never van- 
ishes completely. It was shown with BUU calculations that 
at the balance energy the flow cancellation results from 
a complex transverse momentum dependence and that 
the flow pattern is influenced by EoS and a nn [140]. The 
presently available differential data, measured by FOPI 
down to 90 AMeV [18] suggest that the change of sign of 
V\ is dependent, in addition to transverse momentum, also 
on particle type and rapidity. 

The results on V2 measured at midrapidity are sum- 
marized in fig. 1121 Elliptic flow varies as a function of en- 
ergy from a preferential in-plane, rotational-like [142-144], 
emission (v2 > 0) to an out-of-plane, or 'squeeze-out' 
[37] (i>2 < 0) pattern, with a transition energy of about 
150 A MeV. This transition energy is larger than that 
for the directed flow (see above and the discussion in 
ref. [139]) and was shown to depend on centrality, particle 
type and transverse momentum [28,30]. For higher ener- 
gies, the strength of the collective expansion overcomes 
the rotational-like motion, leading to an increase of out- 
of-plane emission. A maximum is reached at 400A MeV, 
followed by a decrease towards a transition to preferential 
in-plane emission [29, 141]. This behavior is the result of 
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Fig. 12. Elliptic flow parameter V2 at mid-rapidity for colli- 
sions at intermediate impact parameters (about 5.5-7.5 fm) as 
a function of incident energy, in the beam frame. The filled and 
open circles represent the INDRA and FOPI [28] data, respec- 
tively, for Z — 1 particles, the triangles represent the EOS and 
E895 [29] data for protons and the square represents the E877 
data [141] for all charged particles. 



a complex interplay between fireball expansion and spec- 
tator shadowing [28], with the spectators acting as clocks 
of the expansion times. For instance, in the energy range 
400,4 — 1500 A MeV, the passing time of the spectators de- 
creases from 30 to 16 fm/c, implying that overall the ex- 
pansion gets about two times faster in this energy range. 
This interpretation is supported by the observed scaling of 
elliptic flow as a function of transverse momentum scaled 
with beam momentum [28]. We note that the energy de- 
pendence of elliptic flow is similar to that of directed 
flow [1,2,7], with the extra feature of the transition to 
in-plane flow at 4,4 GeV [29] . This high energy transition 
has received particular interest as it is expected to pro- 
vide a sensitive probe of the EOS at high densities [149]. 
At SPS and RHIC energies, the in-plane elliptic flow is 
determined by the pressure gradient-driven expansion of 
the almond shaped isolated fireball [145] and is currently 
under intensive experimental investigation [146-148]. 

The agreement between the corrected INDRA and 
FOPI data is good. The INDRA results have been cor- 
rected using the new method, including the correction 
for the multi-hit losses (see sect. 03 . According to IQMD 
calculations, the reaction plane correction for the lowest 
FOPI energy of 90A MeV appears to be somewhat over- 
estimated. On the other hand, this may partially compen- 
sate for the lack of corrections due to unresolved tracks 
which were not applied for v-i in the FOPI case. Overall, 
the differences between the corrections is small enough, 
so that comparisons of uncorrected data sets are already 
meaningful. A good agreement was found to exist for the 
INDRA [22,33], FOPI [28] and Plastic Ball [37] data in 
the reference frame of the directed flow and without the 
correction for reaction plane resolution [28,33]. 



A remarkable feature of the v-i observable is that it 
allows to show a continuous evolution over a region cover- 
ing completely different reaction mechanisms, from those 
dominated by the mean field near the deep inelastic do- 
main, and the multifragmentation in the Fermi energy 
domain towards the participant-spectator regime at rel- 
ativistic energies. 

7 Correlation between stopping and flow 

Information on stopping and flow in heavy ion collisions 
represents part of the input to theoretical efforts to de- 
duce constraints on the EoS. Remembering that the EoS 
is a relation between pressure and density, it is intuitively 
understandable that these two heavy ion observables are 
related to the EoS: flow is generated by pressure gradients 
established in compressed matter, while the achieved den- 
sity is connected to the degree of stopping. Recently, it was 
observed [7] that a strong correlation exists between the 
stopping, measured in central collisions and the directed 
flow measured at impact parameters where it is maximal 
(see fig. EJ). The relevant data are shown in fig. in the 
upper left panel. Plotted against each other are two dimen- 
sionless global event observables characterizing stopping, 

vartl, and global scaled directed flow, p^Ji r , both defined 
earlier. 

The data points correspond to 21 system-energies with 
varying system size (from Ca+Ca to Au+Au) and energy 
(from 150,4 to 1930,4 MeV). The straight correlation line 
represents a linear least-squares fit to the data and is re- 
peated in the other panels. These other panels show the lo- 
cation along the correlation line of theoretical simulations 
using the IQMD code for Au + Au at 400 A, 1000 A and 
1500,4 MeV as indicated. The points are marked HM and 
SM, respectively, for a stiff (incompressibility K — 380 
MeV) and a soft (K = 200 MeV) EoS. The M in HM and 
SM stands for the momentum dependence of the nn in- 
teraction. IQMD incorporates a phenomenological Ansatz 
fitted to experimental data on the real part of the nu- 
cleoli optical potential. The relevant experimental points 
are given together with their estimated systematic errors 
(these errors were omitted for clarity in the upper left 
panel, but are of comparable magnitude for all the data). 
The main purpose is to show the sensitivity of these com- 
bined observables to variations of the zero-temperature 
EoS as compared to the uncertainty of the data. The EoS, 
that are purely technical, are shown in fig. 1141 The trajec- 
tory of the simulation when changing K seems to follow 
the correlation line, the distance between SM and HM is 
larger at 1000,4 MeV than at 400,4 MeV (i. e. the sensitiv- 
ity is increased at the higher energy), but then does not 
seem to further increase at the highest energy, possibly 
due to the increase of transparency suggested by fig. [7] 
Data measured at energies below 150,4 MeV do not con- 
tinue the same linear correlation, an interesting topic that 
deserves further studies. 

In our exploratory IQMD simulations [97,150] we have 
not tried to be realistic with regard to in-medium modi- 
fications of the nucleon nucleon cross sections a nn , using 
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Fig. 13. Upper left panel: Correlation between the maximal 
directed sideflow and the degree of stopping, after [7]. The line 
is a linear least squares fit to the data, which extend from 0.15 A 
to 1.93 A GeV. The correlation line is repeated in the other pan- 
els which show results of simulations for Au+Au at three in- 
cident energies using two different equations of state, SM and 
HM, together with the experimental points. The short seg- 
ment passing through the SM point in the lower right panel 
shows an estimate of the trajectory using the SM EoS and 
modifying the in medium cross sections in a way that is com- 
patible with [102]. 
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Fig. 14. Equations of state (relative to the ground state) used 
in the calculation. 

instead the vacuum values standardly implemented in the 
code [97]. An estimate of the trajectory in the flow ver- 
sus stopping plot when the EoS is kept constant, but the 
a nn are decreased is shown in the right hand lower panel 
of fig. 1131 For this estimate we rely on the more sophisti- 
cated calculations of ref. [102] which show that a switch to 
more realistic smaller a nn decreases the stopping by about 



20% and the scaled sideflow by 6 — 7%, i.e. a nn acts more 
strongly, relatively speaking, on the stopping than on the 
scaled flow, as expected. The a nn modification trajectory 
crosses the correlation line because it has a different, flat- 
ter, slope than the EoS modification trajectory joining the 
SM to the HM point (which is not plotted) which hap- 
pens to have a slope very similar to that of the experimen- 
tal correlation line. Generally speaking, one can say that 
an underestimation of the apparent transparency will lead 
to an underestimation of the stiffness of the EoS. Never- 
theless, the procedure just outlined suggests that an EoS 
closer to SM than HM would seem to be more appropri- 
ate to describe the data. The same conclusion was reached 
from the comparison of the FOPI data on directed flow, 
including its pr dependence, to IQMD calculations [19] 
and from the comparison of the out-of-plane expansion to 
BUU calculations [32]. 

Despite this encouraging result we would like to stress 
at this time that it would be premature to draw firm con- 
clusions from one particular transport code and it is be- 
yond the scope of this experimental contribution to the 
subject to conclusively settle the question of the EoS. Be- 
sides trying to predict correctly the global observables just 
shown, probably a good strategy to start with, the simula- 
tions must then proceed to reproduce the more differential 
data such as the variations of the stopping and of the vari- 
ous flow components with the particle type, as shown here 
in figs. IFHand llll respectively. Another important physics 
quantity one would like to have under theoretical control, 
in order to be convincing on the conclusion side, is the 
created entropy. Although this is not a direct observable, 
the entropy at freeze-out is strongly constrained by the 
degree of clusterization (of which we showed an exam- 
ple in fig. an d the degree of pionization. An idea of the 
freeze-out volume can be obtained from two-particle corre- 
lations [151], or even multi-particle correlations [152,153]. 
All this is a rather challenging task. We refer to the work 
of Danielewicz, Lacey and Lynch [5] for a summary of 
the situation obtained a few years ago using a subset of 
the then available heavy ion data reaching up to the AGS 
energies. 



8 Summary and outlook 

We have presented a systematics of directed and elliptic 
flow and of stopping for 197 Au + 197 Au reactions in the 
intermediate range of energies from 40 to 1500 MeV per 
nucleon by merging the data from INDRA and FOPI ex- 
periments performed at the SIS synchrotron at GSI. The 
overlap region of the two data sets, 90 to 150 MeV per nu- 
cleon incident energy, has been used to confirm their ac- 
curacy on an absolute scale, and a very satisfactory agree- 
ment has been found. 

Particular emphasis was given to the experimental re- 
construction of the impact parameter and to the correc- 
tions required by the dispersion of the reconstructed az- 
imuthal orientation of the reaction plane. The superiority 
of observables based on transverse energy, either the ratio 
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Erat of transverse to longitudinal energy or the trans- 
verse energy E 1 ^ of light charged particles with Z < 2, 
over multiplicity for the selection of central collisions has 
been demonstrated. A new method, derived by extending 
the Gaussian approximation of the sub-Q- vector distribu- 
tions to the non-isotropic case and by including the effect 
of correlation between the sub-events, has been presented 
and applied to the data at the lower incident energies at 
which the multiplicities are still moderate and the range of 
emitted fragment Z is still large, even in the most central 
collisions. The differences between the standard and the 
new corrections of derived flow parameters are significant 
up to incident energies as high as 150 MeV per nucleon. 

The deduced excitation functions of the v\ and vi ob- 
servables describing directed and elliptic flow exhibit sev- 
eral changes of sign which reflect qualitative changes of the 
underlying dynamics as a function of the bombarding en- 
ergy. The transition from mean-field dominated attractive 
sidewards flow to repulsive dynamics is observed for Z = 1 
and Z = 2 particles at 80 MeV and 60 MeV per nucleon, 
respectively, in mid-central collisions. The transition from 
predominantly in-plane to out-of-plane emissions occurs 
at 150 MeV per nucleon for Z = 1 particles. The sec- 
ond change-of-sign at several GeV per nucleon marks the 
transition to the ultrarelativistic regime. These transition 
points are quite well established and not very sensitive to 
the chosen correction method. The present study shows 
that also the maxima reached by the flow parameters are 
reliable within the typically 5% systematic uncertainties 
due to the corrections and the impact parameter selection. 
Within this margin they may be used to test transport- 
model predictions and their sensitivity to the chosen pa- 
rameterization of the nuclear EoS. 

It has, furthermore, been shown that the significance 
of the comparison can be enhanced by including the exper- 
imentally observed stopping as represented by the ratio of 
the variances of the integrated transverse and longitudinal 
rapidity distributions. This observable can best be deter- 
mined for central collisions at which the directional flow 
vanishes for symmetry reasons whereas the compression in 
the collision zone presumably reaches its maximum. The 
common origin of the observed stopping and flow is ev- 
ident from the strict correlation of the two observables, 
including finite size effects. However, their individual sen- 
sitivity to the magnitude of the nucleon-nucleon cross sec- 
tions and to the flow parameters is different and can be 
used to resolve ambiguities between these two main ingre- 
dients of the models. The sensitivity to parameters of the 
equation of state is shown to increase with bombarding 
energy over the present energy range, and a soft EoS is 
clearly favored by the data. 

Further constraints for the determination of the pa- 
rameters of the equation of state can be obtained by in- 
cluding the detailed dependences of flow on the fragment 
Z, the impact parameter and the accepted ranges of trans- 
verse momentum and rapidity into the comparison with 
theory. These data, for the present reactions, are either 
available already or in preparation. This will have to be 
accompanied by theoretical studies of the still existing sys- 



tematic differences between specific code realizations. The 
importance or necessity of full antisymmetrization at low 
energies or of a covariant treatment at high bombarding 
energies and the role of nucleonic excitations will have to 
be assessed. 

On the experimental side, a gap of missing flow data 
for the Au+Au system exists at energies below 40 MeV 
per nucleon where interesting information on transport 
coefficients as, e.g., shear versus bulk viscosity or thermal 
conductivity may be obtained. The origin of the observed 
negative flow should be confirmed and clarified. At higher 
energies, new information, possibly also on the symmetry 
part of the equation of state, can be expected from new ex- 
periments involving isotopically pure projectiles and tar- 
gets and detector systems permitting mass identification 
at midrapidity. 

The authors would like to thank Y. Leifels for the implementa- 
tion of the IQMD code at GSI, the FOPI and INDRA-ALADIN 
Collaborations for the permission to include partially unpub- 
lished data in this comparative study, and J.-Y. Ollitrault for 
stimulating discussions on flow evaluation and corrections. 
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